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Abstract 

Directional detection of galactic Dark Matter requires 3D reconstruction of low 
energy nuclear recoils tracks. A dedicated acquisition electronics with auto trig- 
gering feature and a real time track reconstruction software have been developed 
within the framework of the MIMAC project of detector. This auto-triggered 
acquisition electronic uses embedded processing to reduce data transfer to its 
useful part only, i.e. decoded coordinates of hit tracks and corresponding energy 
measurements. An acquisition software with on-line monitoring and 3D track 
reconstruction is also presented. 



1. Introduction 

Directional detection is a promising search strategy of galactic Dark Matter. 
The idea is to take advantage on the rotation of the Solar system around the 
galactic center to show a direction dependance of WIMP signals, which should 
then be clearly discriminated from background ones P, 0, 0, 0, [a] ■ Recently, a 
statistical map- based analysis has been developed @ , showing for the first time 
the possibility to extract from data samples of forthcoming directional detec- 
tors, both the main direction of the incoming events, thus proving the galactic 
origin of the signal, and the number of WIMP events contained in the map. 
Several directional detectors are being developed and/or operated : MIMAC 
DRIFT [H, NEWAGE @, DM-TPC [10]. A detailed overview of the status of 
experimental efforts devoted to directional dark matter detection may be found 
in [llj. A common issue amongst these detectors is the fact that directional 
detection requires track reconstruction of recoiling nuclei down to low energies 
(a few keV). This can be achieved with low pressure gaseous detectors [12| and 
several gases have been suggested : CF4, 3 He + C4H10 or CS2. Ideally, recoiling 
tracks should be 3D reconstructed [3| with sense recognition (T3I.JH . 
The MIMAC (Micro tpc MAtrix of Chambers) collaboration 0] is planning 
to build a multi-target detector, composed of a matrix of gaseous micro-TPC 
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detectors. A pixelized bulk Micromegas 15[ is used in order to perform a 3D re- 
construction of few mm tracks. It is segmented in 350 pm pixels associated to a 
dedicated ASIC previously described [Iff . Indeed, the coordinates in the anode 
plane (x and y coordinates) are reconstructed by collecting primary electrons 
produced in the drift region with an electric field (<1 kV/cm) and amplified in 
the avalanche region (20 — » 80 kV/cm). The track of the recoil is thus projected 
on the anode, providing 2D information. As stated above, 3D track reconstruc- 
tion is needed for directional detection of dark matter. This is achieved by 
sampling the anode signal every 25 ns. Knowing the electron drift velocity, in- 
formation on the third coordinate is retrieved. 



To demonstrate the relevance of the concept, a specific acquisition electronic 
has been developed in order to equip a prototype detector featuring an anode 
of 3.36 x 3.36 cm 2 where 2 x 96 strips are monitored with these ASICs. This 
auto-triggered acquisition electronic uses embedded processing to reduce data 
transfer to its useful part only, i.e. decoded coordinates of hit tracks and corre- 
sponding energy measurements. To be fully exploited, an acquisition software 
with on-line monitoring and track reconstruction has been written. 

This paper is organised as follows : section 2 presents the hardware de- 
sign, emphasizing on the front end digital interface (FPGA) implemented on 
the acquisition board. Then the acquisition software and the 3D reconstruc- 
tion strategy are presented in section 3. Eventually, experimental results are 
presented in section 4. 



2. Hardware design and architecture 

2.1. Overview of the MIMAC acquisition board 

As shown on figure [1] the acquisition board is composed of 12 front end 
ASICs, 6 for each coordinate (X and Y). The data processing was split in 3 
FPGA, one for each side and one for merging the data. The LVDS serial links 
are directly connected to the XY FPGA, which de-serializes the ASIC com- 
parators data. Then the XY FPGA processes each ASIC position data for 
suppressing the zeros, decoding the coordinates and assembling them with the 
data coming from the neighbouring ASICs. At the end of the processing, the 
data in the first level FPGA are time sorted and placed in a data buffer (called 
position buffer). 

The two analogue outputs of each ASIC go through peak detectors in order to 
find the maximum amplitude of the shaper signals. Each high gain peak de- 
tector is equipped with a comparator whose threshold is tuned by a DAC. The 
ADC readout is automatically performed by the XY FPGA for each group of 
16 strips when some track has been detected and if the amplitude of the signal 
went above the threshold. As for the position data, the energy data are time 
sorted and placed in a data buffer. 

At this stage, for both channels (position and energy), 2 data sources are avail- 
able: one for the X coordinate and an other for the Y coordinate. The merger 
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Figure 1: Block diagram of the MIMAC acquisition board 



si FPGA is in charge of concatenating them in a time sorted way and make them 

82 available for data readout. An Ethernet ready micro-controller configured as a 

83 TCP socket server is used as a communication link with the acquisition com- 

84 puter. It is at the same time used as the slow control interface for: 

85 • remote configuration of the FPGA (this allows easy upgrades of the firmwares) 

86 • tuning the energy measurement thresholds 

87 • tuning the strips discriminators thresholds 

as In the following, we describe the various components of the acquisition board. 

89 2.2. Front end ASICs 

90 As a preliminary building block for the MIMAC (Micro tpc MAtrix of Cham- 

91 bers) framework a front ASIC was developed [3] . These chips were designed to 

92 monitor 16 strips of pixels with low noise charge pre- amplifiers and to provide 

93 in real time their time over threshold. 

94 The comparator outputs are coupled to a data serializer with a compression 

95 ratio of 8 in order to reduce consumption and connectivity. Indeed, for the 

96 first phase of the MIMAC project, following this prototype phase, the goal 

97 is to equip an 20 cm x 20 cm anode featuring 1024 strips). De-serializing 

98 these data provides a picture of the strips of pixels at a 40 MHz rate. Then, 

99 providing the electron drift velocity in the gas mixture the third coordinate of 

100 the track is reconstructed by using consecutive samples of the strips of pixels. 
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The thresholds are individually tuned by 5 bit current DACs. 
For providing a mean to measure the total energy released in the ionisation 
channel by the incident particle, the sixteen pre-amplifier outputs are summed 
and provided in two gains via two shapers. More details may be found in 16 1. 



105 2.3. Energy measurement with peak detectors 

we Each ASIC shaper output is connected to a peak detector (fig. [5]) which is 
107 designed to memorise the maximum amplitude when enabled. Both high and 
los low gain analogical values are digitized by ADCs when the shaper high gain 
109 output goes back below the threshold programmed by a DAC (comparator part 
not shown on fig. [5]). The first stage amplifies and inverts the negative input 
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Figure 2: Block diagram of the peak detector used for the energy measurement. 



no 

in signal to deliver a positive signal to the second stage. As long as this signal 

n2 increases, the feedback diode stays in conduction mode and the capacitor volt- 

n3 age follows the signal. Once the maximum amplitude is reached, i.e. when the 

n4 shaper signal starts to decrease, the feedback diode switches in blocking mode 

us and thus, the peak amplitude is kept in the capacitor. At this stage, the peak 

no detector amplifier toggles to non linear mode. An anti-saturation loop is imple- 

n7 mcntcd in order to maintain the peak detector amplifier differential input to a 

us minimum value. Indeed one of the bipolar amplifier drawbacks is the fact that 

n9 they exhibit a dramatic input leakage current increase when input differential 

120 voltage increases, which would compromise the analogue memorisation func- 

121 tion. The reset signal is used to disable the peak search and also to discharge 

122 the capacitor at the end of the conversion. Bipolar operational amplifiers were 

123 chosen for implementation as it was required to work with low voltage supplies 

124 and thus having rail to rail input amplifiers was mandatory. 
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2,4' Front end digital interface: XY FPGA 

The XY FPGA0 is composed of two parts. The first one, identical for each 
ASIC, deserializes the data and performs the local triggering, ADC readout and 
local data buffering. The second one is in charge of data time sorting and coor- 
dinate decoding. 

The deserializer is designed to allow operation with low performance FPGA by 
using two shifters working on opposite edges at half of the incoming data fre- 
quency (figure [3]). The sampling phase is adjusted by setting the DCM (Digital 
Clock Manager provided in Xilinx FPGAs) to the appropriate value by slow 
control. 

The ASIC management is detailed on figure @] The deserializer first provides 
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Figure 3: Block diagram of the low cost deserializer 



135 

136 data for the local triggering (an OR of the 16 channels), when a trigger is found, 

137 a copy of the free running time counter is inserted in the raw position memory, 

138 meanwhile the strip data are delayed by one clock cycle to allow this. After the 

139 time count is written in the memory, the strip data are pushed in the memory 
mo while the stop condition is not fulfilled. Experimental results show that primary 
wi electron distribution can be non continuous and therefore strips can be untrig- 

142 gered for a short amount of time. The stop condition is therefore no trigger on 

143 the monitored channel for more than n clock cycles. This value can be adjusted 

144 by slow control. 

145 The memory depth is a key issue. It has to be dimensioned in order to cope 
we with long tracks like alpha particle, because readout is performed at a much 
147 slower speed than the memory writing. For instance, in 3 He at 350 mbar the 
us drift speed is 16 /zm/ns, so an alpha particle arriving orthogonally to the anode 



1 so called because the processing is identical on the X and Y side 
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Figure 4: Block diagram of the ASIC management logic 



149 will use ; — x 40 MHz ~ 375 memory cells. Then, each position word in 

16 jtim/ns 

150 memory can potentially be decoded as 16 valid positions. The memory blocks 

151 used in the FPGA are several times larger than this value. 

152 The trigger is also used by the energy measurement channel. When detected, 

153 the ADC readout system saves the trigger date and waits until a peak is de- 

154 tcctcd within a predefined time limit (shaping time). If none is detected, a null 

155 value is inserted in the energy memory in order to keep raw position memory 

156 and raw energy memory aligned and facilitate later reconstruction. The raw 

157 energy memory is smaller, because for a full particle track only one amplitude 

158 will be coded in two gains. Distributed RAM is sufficient for this task. 

159 Once all data are stored in local buffers, the next step is to assemble them with 

160 the data originating from other ASICs. This operation is done by Finite State 

161 Machine (FSM), like the one used for sorting the position data (figure [5]). The 

162 only difference with the one used for sorting the energy data is the position 

163 decoding (scan word state), each bit at one in position word is decoded as a 

164 strip number. Starting from the idle state, the FSM waits for the first ASIC 

165 (or group) to present data. When it is the case, the current date is saved and 

166 the position decoding is performed. Then the FSM looks if an other group 

167 presents position data marked with the same date, if yes the data are decoded 

168 and appended to the output buffer in the same time slot, if not the time slot is 

169 closed and the search for a new time slot continues. At this stage, 2 possibilities 

170 remain: 

171 • no more data are available and the FSM returns to IDLE 
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172 • data are still available and the earliest data has to be found 

173 For that, a scan of the dates is made (seek_next_group) in order to find the 

174 first ASIC in time. This scanning is complexified by the fact that the time 

175 counter has a short span (13 bit which corresponds to ~1.31 ms), therefore it is 

176 performed in 2 stages. First the search is performed from the current date up 

177 to the maximum counter value, and if not successful, the current date is set to 

178 zero and the searching is performed again until the group is found (check roll 

179 over state). 

lso Position data and energy data are stored in separate buffers (see template in 

lsi table [T]). Both tables are organised in order to optimize bit usage and therefore 

182 data are encoded on 16 bit. For a specific time-stamp, all useful position or 

183 energy data are grouped under a unique TAG. The time slice is closed by the 

184 End of TAG label. 



bit 


15 


14 


13 12 11 


10 9 8 7 


6 5 4 


3 2 10 


TAG 





date 


position DATA n 


1 







ASIC # 


strip f/= 


energy DATA n 


1 





ASIC # 


ADC value 


End of TAG 


1 


1 


unused 



Table 1: Organisation of the 16 bit words containing position or energy information. 



185 2.5. FPGA merger 

186 This FPGA is in charge of merging the data from the X and Y side. For 

187 this task, it is working similarly to the sorting FSM of the XY FPGA, the only 

188 difference being the fact that only 2 data producers are managed instead of 

189 6. Once the data are processed, they are stored in large buffers. They can be 
wo accessed directly for readout by the Ethernet micro-controller. Also, being the 

191 direct interface with the micro-controller, it is used as the slow control interface 

192 for: 

193 • Enabling or disabling an ASIC 

194 • Activating the calibration mode. In this mode a constant pattern is sent 

195 by the ASICs and for reducing the data amount, only one out of 65536 

196 triggers is accepted. 

197 • Configuring the energy channel thresholds 

198 2.6. Ethernet interface 

199 This interface is based on a 32 bit micro controller, using a flash memory 

200 and 16 MB of SDRAM. Its purpose is to provide an easy computer interface 

201 and to test the possibility for multi chamber acquisition. 
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202 At power up, a TCP socket server service is started and the acquisition computer 

203 (the client) can open a socket and communicate with the hardware via a custom 

204 protocol. It is used for slow control and acquisition, they are implemented in 2 

205 threads. When in run mode, the acquisition thread is resumed and start polling 

206 the empty flags of the position and energy buffers, once a non empty flag is 

207 found, the corresponding fifo level is read. Then the appropriate number of 

208 readout is made and the data are transmitted through the socket with a header 

209 specifying the data type and count. 

210 In order to allow flexibility and easy prototyping, the micro-controller firmware 
2n is the only one resident on the electronic board and the 3 FPGA are remotely 

212 reconfigured at each start up via the TCP socket. 

213 3. Acquisition software 

214 The acquisition software was developed in CH — h and coded to allow operation 

215 on Linux and Windows platforms. In order to keep a good separation between 

216 functionalities, it is built with three different parts as described below. 

217 3.1. Software driver 

218 The low level driver is in charge of providing low level functions for con- 

219 trolling the hardware and managing the acquisition. It features a platform 

220 independent non blocking socket client manager, routine for reading FPGA 

221 configuration file and remotely configuring them, methods for interpreting po- 

222 sition and energy threshold files and set DAC values. 

223 For a proper tuning of the serializers sampling phase, a scanning algorithm is 

224 also provided. Basically, it puts the hardware in calibration mode (in order 

225 to significantly decrease the trigger rate), validate one ASIC at a time and try 

226 every possible value of de-phasing (256 values) and record the results of the 

227 reception of a known test pattern. The test pattern is supposed to be received 

228 correctly many thousand times in order to consider that the sampling phase 

229 is correct for the selected ASIC. When the stability range common to the 12 

230 ASICs has been determined, the mean value is computed and transferred to the 

231 XY FPGA. 

232 The data acquisition has to be performed with a blocking call method (returns 

233 only when the required number of bytes are read or when a time out occurs) . 

234 This function reads data frame received from the socket and separates them be- 

235 tween energy and position. They are then pushed in 2 STL (Standard Template 

236 Library) queues. The readout is performed this way, because the electronic pro- 

237 cessing time from particle interaction to data available in computer varies with 

238 energy and position data path. 

239 3.2. Acquisition thread and event building 

240 The acquisition is performed in a thread loop, with the software driver and 

241 Qt framework [17]. The algorithm of the event building is depicted on figure 

242 [Bl When the thread is awaken, it starts looping on the data readout method 
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243 provided by the software driver. As soon as both position and energy queues 

244 are not empty any more the event building processing part starts. 

245 In order to properly assemble an event, the position queue is processed first in 

246 order to find out when it starts and when it stops. All triggered strips with 

247 dates separated by less then the allowed number of empty clock ticks are kept 

248 as member of the same event. But as soon as the date of the currently processed 

249 triggered strips jumps from more than the allowed number of clock ticks, it is 

250 considered to be a member of the next event and the current event position 

251 retrieval is finished. At this stage the event start date and duration are known 

252 and the corresponding energy data may be found in the energy queue. In other 

253 word, the continuous triggering of the anode defines an event. 

254 

255 3.3. Individual strip threshold calibration 

256 In order to operate correctly a MIMAC prototype equipped with the ASIC 

257 described, each strip threshold has to be properly tuned. This task would be 

258 cumbersome if not automated, especially when the aim is to monitor 1024 strips 

259 per chamber. Therefore an automated calibration method was developed and 

260 tested with the prototype. 

261 The procedure is to shut down all high voltage applied to the chamber in order 

262 to have no real events and to try to find the best DAC setting that will put the 

263 thresholds just above noise. This is done by running short acquisitions (<~ 15 s) 

264 with each settings until the good one is found. 

265 First the calibration algorithm will start with all thresholds set at their maxi- 

266 mum values and try to decrease them as long as each individual strip doesn't 

267 trigger (all strip thresholds are tuned at the same time). When a strip has 

268 reached its lowest threshold, it is marked as such. But since there is some 

269 crosstalk on the board, a second phase is necessary. It has been observed that, 

270 when other thresholds lower, some noise is generated on the acquisition board 

271 and picked again by the strips whose threshold were marked as tuned. So in the 

272 second step, the threshold is increased, one digit at a time, until the triggering 

273 is stopped. The chamber is declared calibrated when all strip thresholds are 

274 marked as having reached their minimum value and no more trigger is counted 

275 in the acquisition. 

276 3.4- Graphical interface and event display 

277 As stated above, the readout is performed in 2D with strips of pixels, which 

278 means that an event is considered as valid if, for a given time slot, at least one 

279 strip of pixels in each dimension (X and Y) is fired. When such a coincidence 

280 happens, the position of the fired pixels can be evaluated as the intersection of 

281 the fired X strips and Y strips. Each time slot gives access to a 2D slice of the 

282 track, and each slice following each other will lead to a 3D track. 

283 Once an event is built by the 3D reconstruction algorithm, it has to be 

284 displayed with a visualization software both in online and offline modes. This 

285 real-time visualisation software allows to monitor several parameters of the de- 

286 tector during a run, such as the energy deposit, the position of the interaction 
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287 and the length of the track. Real-time access to these information offers the 

288 possibility to adjust the detector parameters, such as the gas pressure, anode 

289 and micro-mesh voltages or the electronic thresholds in order to reach optimized 

290 working conditions during an acquisition run. Furthermore, the event display 

291 software enables the testing and debugging of the reconstruction software and 

292 the analysis strategy in an offline mode. 



293 
294 



307 



310 
311 



To meet these specifications, a visualization tool mainly based on the Qt 17 1 
and ROOT [l8| framework was developed. Qt offers a well documented user 

295 interface framework implemented in CH — h A third interface layer nammed "Qt 

296 layer" [lj| was used to build the application software. 

297 The user interface of the online software is segmented in several tabs: 

298 • one dedicated only to acquisition process: automated calibration, start/stop 

299 runs 

300 • another tab displays the energy spectrum (both in low and high gain) in 

301 real time 

302 • Event display, used also in offline mode, provides projection of the 3D 

303 track 

304 As shown on figure[71 the main canvas of the event display panel is segmented 

305 in 4 panels: 

306 • the top left panel represents a 2D view of the event, as seen by the anode 



(XY projection) 



308 • the top right panel shows the evolution of the number of fired strips versus 

309 time, for the X side (solid line curve) and the Y side (dashed line curve) 



the bottom pads represent projections of the track on the XZ plane (left 
panel) and YZ plane (right panel) 



312 Event by event display of the projected 3D track (XY, XZ and YZ his- 

313 tograms) is provided by the visualisation software. Whether in online or offline 

314 mode, the list of event to display can be adjusted through several criterion like: 

315 number of images (equivalent to length along Z axis), energy deposit or even a 

316 combination of those. 



317 4. Experimental results of the MIMAC Data AcQuisition 

318 As required by the specifications, the MIMAC DAQ is able to reconstruct 

319 3D tracks of nuclear recoils and electrons between a few keV and hundreds of 

320 keV. 

321 First experimental results were obtained with a 5.9 keV X-ray source ( 55 Fe) 

322 producing 5.9 keV electrons in a 350 mbar 4 He + 5% C4H10 mixture used as 

323 detection medium. These electron tracks (figure [9]) have a specific length over 

324 energy ratio and will be typical background events in the final MIMAC dectector 

325 designed to detect dark matter in an underground laboratory. 
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Ultimately, the goal of MIMAC is to reconstruct recoil tracks of nuclei in 
3D. In order to test the MIMAC prototype associated to the dedicated ac- 
quisition electronics with nuclei recoils, the detector was placed in front of a 
mono-energetic neutron beam at the Amande facility 20]. Depending on the 
target gas used, recoils of hydrogen ions, helium ions or event fluorine ions were 
reconstructed [21]. 

For instance, on figure [8j a 100 keV proton recoil obtained in a 350 mbar 
4 Hc + 5% C4H10 mixture is displayed in three dimensions. Left panel presents 
the raw 3D track, ie images of the anode in each z slice. Right panel presents 
the barycentre of each z slice. This shows the possibility to achieve 3D track 
reconstruction of recoiling nuclei in low pressure gaseous TPC. Figure[5]presents 
the same views of 5.9 keV electrons, highlighting the possibility to reach low 
energies with this technique. This last event is of crucial interest as it features a 
typical background event for Dark Matter search. A comprehensive study of an- 
gular resolutions obtained with this 3D reconstruction method will be presented 
in a forthcoming paper 22]. This dedicated acquisition electronics associated 
with 3D reconstruction software offers a major breakthrough towards 3D recon- 
struction of low energy tracks, thus opening great opportunities for directional 
Dark Matter search 12 ll. 



345 5. Conclusion 

346 A dedicated acquisition electronics with auto triggering feature and a real 

347 time track reconstruction software have been developed within the framework 

348 of the MIMAC project of detector. It has been shown by experimental results, 

349 that the MIMAC ASIC [ll], the prototype acquisition board (fig. [TUJ) and soft- 

350 ware coupled together offer great possibilities for 3D track reconstruction and 

351 ultimately for directional dark matter detection. 
352 
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Figure 7: Snapshot of the event display tab of the MIM AC real-time visualization software. 
Top left panel represents a 2D view of the event, as seen by the anode (XY projection). Top 
right panel shows the evolution of the number of fired strips versus time, for the X side (solid 
line curve) and the Y side (dashed line curve). Bottom pads represent projections of the track 
on the XZ plane (left panel) and YZ plane (right panel) 



14 




Figure 8: 100 kcV H recoil in a 350 mbar 4 He + 5% C4H10 mixture: left panel represents the 
raw 3D track while the right panel shows barycentres of each time slice 




Figure 9: 5.9 keV electron recoil in a 350 mbar 4 He + 5% C4H10 mixture (barycentre on the 
right panel) 
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Figure 10: Picture of the 260 mm X 246 mm MIMAC acquisition board. Ethernet interface 
board and connectors are missing on this view. X side of the dectector is connected on the 
top connector and Y side on the bottom one. The 2x6 ASIC are fed by the connectors. In 
the middle the 3 FPGA perform the processing. 
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